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Abstract 



This paper is devoted to the study of the Cauchy problem for the Boussinesq system 

with partial viscosity in dimension N > 3. First we prove a global existence result for 

data in Lorentz spaces satisfying a smallness condition which is at the scaling of the 

^2 ■ equations. Second, we get a uniqueness result in Besov spaces with negative indices of 

regularity (despite the fact that there is no smoothing effect on the temperature). The 

proof relies on a priori estimates with loss of regularity for the nonstationary Stokes 

system with convection. As a corollary, we obtain a global existence and uniqueness 

result for small data in Lorentz spaces, 
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^ ! 1 Introduction and main results 

00 . 

The present paper is devoted to the mathematical study of the so-called Boussinesq system 
with partial viscosity: 



d t 6 + u-V9 = 

dtu + u-Vu- vAu + VTl = 8e N (1) 
x divu = 0. 

Above, 9 = 9(t,x) and n = IL(t,x) are real valued functions, and u = u(t,x) is a time 
dependent vector field. We denote by e^ the vertical unit vector of M. N . It is assumed that 
the space variable x belongs to R . We supplement the system with Cauchy conditions 
(9o, uq) at time t = and address the question of solvability for t > 0. 

Boussinesq system arises in simplified models for geophysics in which case u stands 
for the velocity field and the forcing term 6 e^ is proportional either to the temperature, 
or to the salinity or to the density (see [17]). Here, we shall call 9 the temperature. 

Remark that the standard incompressible Navier-Stokes equations arise as a particular 
case of ([!]) (just take 9 = 0). Hence, it is tempting to study whether the classical theory 
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for the Navier-Stokes equations may be extended to those more general fluids. As far 
as one is concerned with global results, the main difficulty that we have to face is that 
9 is transported by the flow of u. Hence it has no time decay nor gain of smoothness 
whatsoever and the standard approach for solving Navier-Stokes equations with a (given) 
source term is bound to fail. Nevertheless, in a recent paper (see [H]), we have stated 
that system (pQ) has global finite energy weak solutions for any data (9q,uq) in L? and 
that uniqueness holds true in dimension two. That result may be compared with Leray's 
theorem [16j for the Navier-Stokes equations. 

In [11], we have also obtained a global existence result (for small data) in the spirit of 
Fujita and Kato's theorem (see [12] and p3]) for the Navier-Stokes equations: 

Theorem 1 Let N > 3. Let 9 £ Lis (R N )nB° NA (R N ) and u € L N '°°(B. N )nB pA v (R N ) 
for some p £ [N, oo] . There exists a positive c depending only on N and such that if 

\\Uo\\ L N,aa + 1 1 1 6*0 1 1 — CV 

then Boussinesq system ([1]) admits a unique global solution 

(M,VII) eC(R+;B° Nil ) x (c(K + ; J B| 1 _1 )n^L(K+; J B| 1 +1 )) JV x (^(R+j ^JT))'" . 

Besides, from dimension N = 4 on, space Lt(M. n ) may be replaced by Lt ,00 (M w ). 

The reader is referred to section [2] for the definition of Lorentz spaces L p,oc and Besov 
spaces Bp r which have been used above. Remind that we have the (continuous) inclusions 

B° N>1 C L N C L N '°°. 

Theorem Q] may be interpreted in terms of scaling of Boussinesq system. Indeed, ([1]) is 
obviously invariant under the transform 

u(t, x) i — ► Xu(X 2 t, Xx) and 9(t, x) i — > X 3 6(X 2 t, Ax) for all A > 0. 

Hence data (^o ; ^o) belong to a functional space E which is at the scaling of the system 
if and only if the norm of E is invariant by 

Uq{x) i — > Xu Q (Xx) and o (x) i — > X 3 9 (Xx). (2) 

In dimension > 3, the spaces if x (L N >°°) N and if' 00 x (L N >°°) N satisfy ©. Hence, 
Theorem Q] is a global well-posedness result for suitably smooth data satisfying a scaling 
invariant smallness condition. In fact, prescribing additional regularity in homogeneous 
Besov spaces ensures that the velocity belongs to the set L z 1 oc (IR_|_; Lip) of locally integrable 
functions in t with values in the set of Lipschitz vector fields. This property was needed for 
uniqueness. Note in passing that the additional condition required for the initial velocity 
still satisfies ([2]). 
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The present paper aims at weakening as much as possible the Besov space assumption. 
Ideally, we would like to consider general (small) data (60, uq) in the scaling invariant 
space lf-» x (L N '°°) N . 

For proving existence, the strategy is the following. We smooth out the data so that 
Theorem [T] provides a sequence (9 n ,u n ) of global solutions which is bounded in 

L°°(R + ; L^"' 00 x (L N '°°) N ). 

Obviously, those bounds are insufficient to pass to the limit in the nonlinear terms. How- 
ever, the parabolicity of the second equation of ([I]) provides additional regularity so that 
there is some chance to pass to the limit anyway provided 9q belongs also to L p '°° for 
some large enough p. Those basic considerations will enable us to prove the following 
statement: 

Theorem 2 Let uq be a solenoidal vector field with coefficients in L N '°°. Assume that 6$ 
belongs to 

• L 1 n L p '°° for some p>3/2 if N = 3, 

• L 4 / 3 ' 00 n LP' 00 for some p > 4/3 if N = 4, 

• L N / 3 '°° n L p '°° for some p > N/3 if N > 5. 

There exists a constant c > depending only on N, and such that if 

\\uq\\ L 3,oo + H^oIIl 1 < cv if N = 3, ||uo||iiv,oo + 1 1 6*o H^/s.oc <cv if N > 4 (3) 
then system ([1]) has a global solution (0, it, VII) with u G L°°(R + ; L N,oc ) and 

G L°°(1L|_; L 1 n L p '°°) if N = 3, 9 G L°°(IR + ; n L p,co ) if N > 4. 

Remark: In fact, the heat semi-group supplies some additional regularity properties for u 
(that we shall use in the proof of Theorem [2]) . 

Note that from dimension 5 on, one can take p = N/3 so that the above statement 
is a global existence result in a scaling invariant space for the system. The scaling of the 
system may still be almost achieved in dimension 4. In dimension 3 however, it is very 
unlikely that an existence result may be proved for #o £ L l and no G L N >°° ( se e Remark [4] 
for further explanations). 

As regards uniqueness, let us stress that for general no in L N '°° the function e* A no 
(where (e tA )t>o denotes the heat semi-group) need not be in Lj oc (M. + ; Lip). Therefore, the 
velocity field of the solution constructed in Theorem [2] need not be in Lj oc (M. + ; Lip) either 
which precludes us from proving stability estimates for system (pQ) by mean of standard 
arguments. Indeed, we have to deal with a transport equation associated to a vector field 
which is not in Lj oc (JH + ; Lip). This difficulty has been overcome in [TT] in the framework 
of two-dimensional finite energy solutions. In the present paper, we shall see that similar 
arguments may be used to state uniqueness in dimension N > 3. 
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Theorem 3 Let {9\,ux, VIIi) and (#2>it2, VII2) satisfy (pQ) wit/i f/ie same data. Assume 
that for some p E [1, 27V[ and i = 1,2, we have 

-l+E -1+^ ~ l+— 

^ € L§?(B Pl00 P ) and Uj G L??(5 Pi00 p ) n L^(B Pi00 p ). 

There exists a constant c > depending only on N and on p such that if in addition 

\\ui\\_ . 1+ JV +U~ 1 \\U2\\ .-i+JV <c (4) 

then VIIi) = (# 2 ,u 2 ,Vn 2 ). 

In the above statement, the space L\,{B Pt0 ^ ) is slightly larger than the set L)p(Bp j00 p ) of 

integrable functions over [0, T] with values in the Besov space -B p ,oo p (see section [2]). 

Remark: The limit case p = 2N may be considered if the velocity field belongs to a 
Besov space with third index 1 (see Theorem [5]) . 

■ -1+— 

Finally, putting together the embedding L N,oc «— > B qtOQ q for q > N (see Lemma [2] 
below), the existence and uniqueness theorems, and the further regularity properties for 
the velocity given by the heat semi-group, one ends up with the following global well- 
posedness result: 

Theorem 4 Assume that (9q,uo) satisfies the assumptions of Theorem with p = N. 
Then system ([I]) has a unique solution (6,u, VII) such that u G L°°(M.+ ; L N,oc ), 9 G 
L°°(R+; LT' 00 ) if N > 4 and 9 G L 1 (R+; L ) if N = 3, with moreover 

~ 1+^- 
ueL] oc {B q ^) for all q>N. 

The paper is structured as follows. In the next section, we present a few tools borrowed 
from harmonic and functional analysis. Section [3] is devoted to the proof of existence. The 
study of uniqueness is postponed in section [H 

Acknowledgments: The authors are indebted to the referee for fruitful remarks. 

2 Tools and functional spaces 
2.1 Lorentz spaces 

To start with, let us recall the definition of weak LP spaces (denoted by L p, °°): 

Definition 1 For 1 < p < 00, we denote by L p,oc (M. N ) (or simply L p '°° ) the space of all 
real valued measurable functions over W N such that 

1 

< 00. 



LP- 



sup A 



{xeR N /\f(x)\>\} 
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Remark 1 The space L p '°° may be alternatively defined by real interpolation: 

L^ = (L-,L 1 ) ( i i0o) . 

In other words, a function f belongs to L p '°° if and only if, for all A > 0, one may write 
f = f A + fA f or some functions fA £ L 1 and f A in L°° such that 

\\f A \\ L i KCA 1 ' 1 ^ and \\f A \\ L ~ <CA- l /v. 

The "best constant" C defines a norm which is equivalent to ||/||£p,°o. 

The set C£° of smooth compactly supported functions is not dense in spaces L p '°°. It 
turns out however that is locally dense in L p '°° (despite the fact that L p '°° is not 
separable). More details are given in the following proposition. 

Proposition 1 For all p g]1,oo[ and f G L p,oc , there exists a family (f £ )e>0 °f 
functions and a constant C so that 

sup ||/ £ || L p,oo < C||/||lp,°o and f £ ^f in L 1 + L°° . (5) 

£>0 

Proof: Let ((p e )e>o be a family of mollifiers and (Xe)e>0i a family of cut-off functions 
with values in [0,1], supported in 5(0, 2e _1 ) and equal to 1 in a neighborhood of 

5(0, e- 1 ). 

Let / £ L p '°°. For all e > 0, set f £ := (p E * (xef) an d fix some rj > 0. From the 
above definition and remark, it is obvious that 

Wfehp^ < C\\xef\\LP,o° < C||/||lp,oo for all e > 0, 

and that one can find two functions g € L 1 and h € L°° such that / = g + h and 
< 77/4 (just take A large enough). 

Let us split f e — f as follows: 

fe~ f = {(fe* lXe9) ~ * ^} + \fe * 9 ~ ff} + ★ (x e h) ~ h\ . 

On the one hand, Lebesgue dominated convergence theorem and standard results 
on convolution ensure that the first two terms between curly brackets have L 1 norm 
less than rj/2 for small enough e. On the other hand, the L°° norm of the last term 
is obviously less than rj/2 for any e > 0. 

It is now easy to complete the proof of the proposition. ■ 
2.2 Besov spaces 

In this subsection we define the Littlewood-Paley decomposition and the Besov spaces 
that we are going to work with. The reader is referred to the monographs [5], [15] or [18] 
for a more detailed presentation. 
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To start with, fix a smooth nonnegative radial function x with support in the ball 
{|£| < §}, value 1 over {|£| < |}, and such that r i— > x( re r) be nonincreasing over R + . 
Let </?(£) = x(£/2) — x(£)- We obviously have 

v{2~ q t) = 1 for all £eR N \ {0}. (6) 



We define the spectral localization operators A q and S q (ggZ) by 

A g u := ip(2~ q D)u and S, it := x(2~ 9 £>)u. 

For any tempered distribution it G cS'(R JV ), functions A g it and S^it are analytic with at 
most polynomial growth and u = A q u modulo polynomials (see e.g. |15j). 



We shall often use the following quasi-orthogonality property : 

A k A q u = if \k-q\>2 and A k (S q -iuA q v) = if |&-g|>5. (7) 
Let us now recall the definition of homogeneous Besov spaces. 
Definition 2 Let s G R, (p,r) G [1, cx)] 2 and it G ^'(R^). We set 



i 

, || / X ~ r^rriQU \ 1 1 1~ \ 

it 



i / z A„u rJ 2/ r < oo, and it a s := sup 2 q A ? u 



If s < ^ or s = — and r = 1 i/ten f/te homogeneous Besov space Bp r := ^ ^(R^) 
is defined as the set of those tempered distributions u such that \\u\\^ a < oo. 



• Iff + k<s<f + k + l (or s = f + k + 1 and r = l) for some /c G N then B^ r 

is the set of tempered distributions u so that d a u G Bp 9 ^^ 1 for all multi-index a of 
length k + 1. 

Remark: Let us all also recall in passing that the nonhomogeneous Besov space Bp r is 
the set of tempered distributions u such thalQ 



|u||b s . < oo with IIwIIb* , := ||<5qIi||zj> + 2 9S ||A„it| 



LP 



Remark: For all s G R, the Besov space B\ 2 coincides with the homogeneous Sobolev 
space H s , and, if s G R + \N then B^ ^ coincides with the homogeneous Holder space C s . 

The following Bernstein inequality will be of constant use in the paper. 

Lemma 1 Let k G N, 1 < p\ < p2 < oo and ip G C%°(M N ). There exists a constant C 
depending only on k, N and Supp^ such that 

\\D k iP(2- q D)u\\ L P 2 < C2 q ( k+N ^~^\\%l,{2~ q D)u\\ L v 1 . 



1 note that the definition is somewhat simpler than in the homogeneous framework since low frequencies 
cannot experience divergence. 
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We shall also use the following fundamental properties of Besov spaces. 



Proposition 2 (i) The set B^ r (R ) is a complete subspace of <S'(R ) if and only if 

s < N/p or s = N/p and r = 1. 

(ii) There exists a positive constant c so that 

c~ 1 \\u\\f>s < II Vull ns-i <c||u||p s . (8) 

D p,r "&p,r P,r 



J y Pl VI' 

3 P2,T2 



(Hi) For all 1 < p\ < p2 < oo and 1 < rx < r*2 < oo, we have B^ iri <— > B p 

JV JV JV 

(iv) If p € [1, oo] then 5 p 1 =— > -Bp^oo H // p is /mite then space B^ is an algebra. 

(v) Real interpolation: (B^ r , i?*^)^, = B 6 ^^ 1 9 ^ Sl whenever < 9 < 1 and 1 < 
Pj r, r' < oo. 

We shall often use the fact that Lorentz spaces are embedded in Besov spaces (see the 
proof in 

Lemma 2 For any 1 < p < q < oo, we have 

JV JV 

L P,oo( R jV) ^ ^ P ( R jV)_ 

In order to pass to the limit in the nonlinear terms of System ([I]) , the following compactness 
result in nonhomogeneous Besov spaces will be most useful. 

Lemma 3 For any 1 < p < q < oo, e > and ip £ C£°, the map u i— > ipu is compact 

JV JV £ 

/ram L^QR^) to B q % p (R N ). 

Proof: It is well known that u i— > ^/m is a compact mapping from Bp r to Bp' r whenever 
s' < s (see e.g [18] and the references therein). So it suffices to combine the previous 



lemma with the fact that, owing to N/q — N/p < 0, we have B q q )0 o v ^ B q q j00 



2.3 Product and paraproduct in Besov spaces 

In order to bound the nonlinear terms in system (pQ), it will be useful to know how the 
product of two functions operates in Besov spaces. In fact, optimal results may be achieved 
by taking advantage of (basic) paradifferential calculus, a tool which was introduced by J.- 
M. Bony in [4j. More precisely, the product of two functions / and g may be decomposed 
according to 

fg = f f g + f g f + R(f,g) (9) 
where the paraproduct operator T is defined by the formula 

and the remainder operator, R, by 

R(f, g) := N/A 9 5 with A q := A^i + A q + A^. 

q 



We shall make an extensive use of the following results of continuity for operators T and 
R (see the proof in e.g. [18]): 

Proposition 3 Let 1 < p,pi,P2,r,n,r 2 < oo so that | = i + i and i = i + i. 

Operator T is continuous : 

• /rom L°° x B* to B* /or a// teM, 

• from B~* ri x B* 2r2 to B*~ s /or all t £R and s > 0. 
Operator R is continuous: 

• /rom B^ x to B*+* /or all (s,t) G M 2 swc/i f/iat s + 1 > 0, 

• from B s pun x B~% to B$ >00 if s eR and ± + ± > 1. 

Combining the above continuity results with Bony's decomposition ([9]), we get: 

Proposition 4 Let (p,r) € [l,oo] 2 and (s\,S2) € M 2 . T/ie following inequalities hold 
true: 

• If si + s 2 + JVinf(0, 1 - 2) > 0, si < f and s 2 < f then 



uv 



p , V ^^p j T* ^ 



If s i = — then \\u\\ n (resp. \\v\\^ s 2^ ) has to be replaced with ||it|| jv (resp. 

N 



IV 



|o s 2 ). If S2 = then \\v\\ jv /ias to 6e replaced with \\v\\ n 

a P,r ' V R P R P 



R P R P nL c 



// si + s 2 = 0, si e (-f, f] andp>2 then 



If p > 2 toen 



|HI w < Gllull^x ||u||^ . (11) 

D D . OO 



\uv\\ _n < C\\u\\ _n\\v\\ n . (12) 

B p B v B p nL°° 



// \s\ < y and p > 2 or —y < s < ~ and p < 2 then 



\uv\lfts <C||n||r> s \\v\\ jv . (13) 
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2.4 Estimates for the heat and Stokes equations 

We shall often use the following smoothing property for the heat equation which has been 
stated by J.-Y. Chemin in [6]. 

Proposition 5 Let s G R, 1 < p, r, p\ < oo. Let uq G B* and f G L\f {B Pyr P1 ). Then 
the heat equation 

d t u - vAu = f, tt| t=0 = u 

~ . ~ -s+j- 

admits a unique solution u in L^(B pr ) fl L^ (B P:1 . Pl ) and there exists a constant C 

depending only on the dimension N so that for all t G [0, T] and p > p±, we have: 

HML . s+ | <c(\\u \\ 6 +^~ 1 II/IL Y (14) 

L p t {B P ,r P ) V p ' r L P1 {B P ,r P1 )J 

In the above statement, spaces Lj,(B pr ) are defined along the lines of definition [2] with 

( 15 ) 



^II^'wIIl^Lp) 



Note that by virtue of Minkowski inequality, we have 

IMIz£(%„) < W U ke T (B°, r ) for P ^ T > ( 16 ) 
and the opposite inequality if p < r. 

We shall set Lf oc (R+; B^ r ) = f] T>Q L p T (B s pr ). 

Remark 2 The above results of continuity for the paraproduct, remainder and product 
may be easily carried out to L P r ,(B pr ) spaces. The time exponents just behave according 
to Holder inequality. 

Remark 3 Since the Leray projector V over solenoidal vector fields maps B p r to B pr 
for every s£l and 1 < p, r < oo, Proposition [5| may be extended to the nonstationary 
Stokes system 

dtu — vAu + VII = /, divu = 0, . . 

u\ t=0 = n , 

with divergence free initial data uq G B pr , and source term f in L^p(B p r ). In particular, 
denoting Q = Id — V, we have the following a priori estimates for all p>\ : 

^ IH W(C f ) - + and l|Vn|l 4(^) - c ^h(Bi r y 

3 Existence of solutions with infinite energy 

This section is devoted to the proof of Theorem [2j The principle of the proof is standard: 

1. approximate the data (6o,uq) by a sequence (0Q,UQ) n ^ of smooth solutions; 

2. solve (globally) the Boussinesq system with data (6q,Uq); 

3. resort to compactness arguments to prove the convergence of a subsequence; 

4. pass to the limit in the system. 
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3.1 Global existence: the smooth case 

As a warm up, let us first consider the case p = N which is easier to deal with. For 
notational simplicity, we agree that L"3"'°° stands for L if N = 3. 



Let (0 O , uq) be in (if' 00 n L^' 00 ) x (L N >°°)" with dWu = 0. According to (a slight 
generalization of) Proposition [Q one can find a sequence (#o,Uo)«eN of functions 
which tends to (9q, uq) in the sense of the distributions and such that in addition 

ll^o ll L f,oo < Cll^oll^^.oo, II^oIIl^ 00 < CII^oIIljv.oc, HmoIIl^ 00 < C\\ u o\\l n .°°- 

Note that projecting on the set of solenoidal vector fields is needed to ensure that divtiQ = 
0. This operation is harmless since V maps L N,oc in L N,oc . 

According to [11] Theorem 1.5, we deduce that there exists a positive constant c so 
that if condition ([3]) is satisfied then system ([1]) admits a unique global solution (6 n ,u n ) 
in C(M+;i?jY 1 ) with in addition 



u 



in,oo<C{v 1 ||0o|| n ^ + \\uo\\ l n,oc) for all t > 0. 



As 9q G L 3 >°° n L '°° and divu n = one can assert (see e.g Proposition 4.6 of [IT]) that 

||0 n (*) || n i00 < C\\e \\ n >00 and \\e n {t)\\ LN ,ao <C\\e \\ L N, x forall t > 0. 

Note that the above estimates are not sufficient to pass to the limit in the nonlinear terms 
of system ([TJ). So we shall take advantage of the smoothing properties given by the Stokes 
operator in the equation for the velocity. Let q be in ]N, oo[. Putting together Remark [3] 

and the embedding L N,oc B q q :00 , we discover that 

^|K"|L .£+1 <C(\\u Q \\ L N >a o + \\li n ®U n \\„ t N +11^111,1(1,^,00)]. (18) 

Using Bony's decomposition ([9]) followed by Proposition [3l Remark [2] and Lemma [2j we 
can write that 

N 

< C \\ u roo/R" 1 \\\ u Kj-i > 

<_ C|| U 71 1 1 £00 (£,JV,oo) || It 71 1 1 _ 



Plugging this latter inequality in (|18p . and assuming that c is suitably small, we thus get 
^IKIL .js+i < c(ll moIIl^.oo + ^||^o||ljv>o°) • (19) 



Now, compactness arguments will enable us to pass to the limit. Indeed, we have proved 
that sequence (6 n ,u n ) n< =n is bounded in 



L°°(1R + ; LT -°° n L N,co ) x f L N,oc ) n LL(M+; ) ) for all q > N. 
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So it is easy to show that sequence (<9t# n ) ng N is bounded in the set of space derivatives 



N _ „, -3+ 



N 



of functions of L°°(R + ;L 2 >°°), which is embedded in L°°(IR + ; B qt00 q ). Now, according 



3+ 



AT 



to Lemma El for all tp G C£° the map / 1— * ipf is compact from L >°° to Bq )00 q . So one 
may conclude by combining Ascoli theorem, Cantor diagonal process and interpolation 
that there exists some function 9 G L°°(IR + ; L~ n L N,oc ) such that, up to extraction, 

iP6 n ^ipe in Lf^.(R+; Bg t00 9 ) for all ^C c °°, ee]0,l[ and q > N. 

N 

Similar arguments show that (u n ® -u n ) ne N is bounded in the space L~'°°) and 

N 

that (0 n ) nS N is bounded in L°°(R + ; L~ ,0 °) so that, using embeddings, we conclude that 

((9t^ n )neN is bounded in L°°(M + ; i? gj00 9 ). Repeating the above compactness argument, 

~ ■ ^-+1 

we get some distribution u G L°°(IR + ; L N '°°) n L\ oc {W + ; B q q oo ) so that, up to extraction, 

-l+£-e 

iPu n ^ipu in L% C (R + ; B q ,oo q ) for all ip G C^°, £G]0,1[ and g > iV. 

~ - i+£ 

Interpolating with the uniform bounds in Lf oc (M.+ ; B q:00 q ), we discover that convergence 

for (ipu n ) ne ^ also holds in every space L[ oc (R + ; B q q i00 r ) with r > 1 and q > N, which 
suffices to pass to the limit in all the nonlinear terms. So (9, u) is a (weak) solution to 
system (JT]). ■ 

3.2 Global existence: the general case 

Let us now prove the existence part of Theorem [2] in the general case. Let n G L N >°° and 

N 

9q G L"3"'°° n L p,oc for some p satisfying the conditions of Theorem [2j Assume that the 
smallness condition fl3J) is satisfied. 

As before, we solve system (pQ) with smoothed out data and obtain a solution (9 n ,u n ) 
in C(M + ;i?^ ri ) satisfying 

\\u n (t)\\ L N,oo < c(^ _1 ||6»o|| L ^ i00 + HuoHiiV.oo), 

\\e n (t)\\ L ^ < \M L ^ and \\9 n (t)\\ LP ^<\\9 \\ LP ^. 



Note that u n satisfies 



u n (t) = e tuA u% 



I e^- T > A V div (u n ® u n ) dr + I e (*—> A p(fl™ e N ) dr . 
Jo Jo 



By virtue of the embedding L N,oc <— > Bql iOC for q\ > N and of Proposition [5l we get 
~, . f +1 

G L loc (K + ; B q ^oo ) uniformly with respect to n for gi > N. (20) 
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Next, we notice that (u 11 ® u n ) n ^ is bounded in L°°(R + ; L 2 >°°) so that, combining the 

N ■ ^ 

embedding L~ ,co Bq^oo for 52 > iV/2 and Proposition [5l 

«2 e^ c (K+;Sg 9 2 2 j00 ) uniformly with respect to n for g 2 > iV/2. (21) 

AT _ JV 

Finally, (6> n ) neN is bounded in L°°(M + ; and L p '°° <^-» J B g 9 |, 00 p for g 3 > p hence 

jv _ jv_i_2 

«3 eLi£ c (R + ;Bqloo P ) uniformly with respect to n for g 3 > p. (22) 

Therefore, compactness arguments similar to those which have been used in the previous 
section enable us to show that, up to extraction, (8 n ,u n ) n ^ tends in the sense of distri- 

N 

butions to some (6,u) such that 6 G L°°(R + ; L~ ,co PiL p,ao ) and, for all q x > N, q 2 > N/2 
and q3 > p, 

u G L°°(IR + ; L N n n (lL(»+; <,oo ) + (M+; B£ >00 + i^ 3 ,co P )) • 

By interpolation with the uniform bounds stated above, we deduce that convergence holds 
true for all e G]0, 1[, q > p, q\ > N, q 2 > N/2 and q% > p, 

— — — — e 

• locally in Lf£ c (R + ; B q q )OQ p ) for the temperature, 

• locally in L r (R + ;B^ )OQ r ) + L% C [R+; Bg t00 + B^ ) for the velocity. 

Taking advantage of continuity properties for the paraproduct and the remainder, it is 
then possible to pass to the limit in u"# ra , u 2 9 n and u^8 n whenever there exist some 
exponents q > p, q\ > N, q 2 > N/2 and q%> p satisfying 

K_n + n + 1 jVmax(o,± + ^--lY (23) 
q p qi \ 7 q qi J ' v 1 

K_n + n_ 1 Nmax ( 1 + X_ 1 ) (24) 

q P ?2 V 9 92 / 

K_ n + n _ N +2 jVmax(0,± + ^--lY (25) 

It is clear that the first condition is satisfied for q close enough to p. Next, second condition 
is verified if q 2 is close enough to N/2 provided p > N/(N — 1). Finally, one can find 
some (73 so that the third condition be fulfilled if and only if p > j^t^. 

Note that in dimension three, this implies that p > 3/2, and that in dimension four, 
we must have p > 4/3. From dimension five on, one may find some q, q±, q 2 and q^ such 
that conditions (I23p . (I24D and (J25]) are satisfied with p = N/3. That no further condition 
is needed to pass to the limit in u n ®u n is left to the reader. As a matter of fact, because 

m"0/ = u"0 « + u% + u%), 

it suffices to put together the fact that (u n ) n& fq is bounded in L^ C (M + ; L N,oc ) and the 
properties which have been stated above for sequences (u™) ng N, (u 2 ) ne n and (u^) ne ^. ■ 
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Remark 4 In dimension three, one may prove global existence under the weaker condition 
that p > 6/5. This may be achieved by combining the bounds satisfied by (u") n gN> (w^neN 
and (u^neN with a bootstrap argument. 

It is not clear however that one can take p < 6/5. Indeed, having 9 n in L°°(IR + ; L P) °°) 

implies that is no better than L^ C (M + ; Bq^ >00 p ) for all > p. Hence p > 6/5 is 
needed to pass to the limit in div (6 n u^). 



4 Uniqueness 

Let us first give the heuristics of the proof of Theorem [3j 

Consider two solutions (9i,u%,VILi) and (02, U2, VII2) of system ([1]) corresponding to 
the same initial data. Assume that for some 1 < p < 2N, 

-1+— / — 1+— ~ 1+— \ N 
(Pi,iH) € Lf(B Pt00 p ) x (Lf(B p>00 p ) n L X T {B P ,J )J for t = 1,2. 

The system satisfied by the difference (59, 5u, V&I) between the two solutions reads 
f + div (uid3) = - div (9 2 6u), 

\ dtdu + div (ui ® 5u) — vASu + V<fl = — div (5u (8 W2) + $ e-N- 

Note that the right-hand side of the first equation (which is a transport equation asso- 
ciated to the vector field u\) is (at least) one derivative less regular than 62. Because 
no smoothing property may be expected for such an equation, this obliges us to perform 

estimates in L^ } (B PtOC p ) for 59 rather than in L^(B Pt00 p ). Now, due to the coupling 
between the equations for 5u and 59, this loss of one derivative also occurs in the esti- 
mates for 5u. This yields the constraint p < 2N when bounding the quadratic terms (see 
Proposition S]). 

The second difficulty that we have to face is much more serious: as the space L^(i? Pi00 p ) 
fails to be embedded in L^(Lip), the vector field u\ is not in L^(Lip). Therefore, the 
initial regularity of 59 need not be preserved during the evolution. It turns out however 

that B Pt oJ is embedded in the set Loglip of Log-Lipschitz functions so that one may 
resort to arguments similar to those used by H. Bahouri and J.-Y. Chemin in [3] to prove 
estimates with (small) loss of regularity. Of course, we will have to cope with the fact that, 

due to the "tilde", the space L^(i? P)00 p ) is not quite embedded in L^(Loglip). Overcoming 
this ultimate difficulty is the purpose of the next subsection. 



4.1 A priori estimates with loss of regularity 

This section is devoted to proving a priori estimates with loss of regularity for transport- 
diffusion equations of the type 

(d t p + div (pu) - uAp = f ^ 
I p\ t =o = po 
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with u a given solenoidal vector field with coefficients in L^(£? P)00 P ) or, more generally, 
for the nonstationary Stokes system 

d t v + div (u ® v) - vAv + VII = /, 
div u = 0. 

Let us first state a result for equation (|26|) which is a trifle easier to deal with. 

Proposition 6 Ze£ 1 < p < oo and s G ] —1 — min(y , =7) , 1 + [. Lei p be a solution to 
the transport- diffusion equation (|26|) . There exists some Nq E N depending only on the 
choice of the Littlewood-Paley decomposition, a universal constant Co and two constants 
c and C depending only on s, p and N so that if 

llVdl K <c (28) 

rlfD P \ 

T v p, 00 y 

/ iV 

i/ien := C \~] 2 q ( 1+ T^ / ||A 9 'm||/,p dr with A q i = Yl\ a \<N Ag'+a satisfies 

Jo 

e q (T)-e q/ (T)<±(l + s + min[f,f))(q-q') for all q > q> 
and the following a priori estimate is satisfied for all t G [0, T] : 

sup 2 qs - £ ^\\A q p{T)\\ L v +^sup / 2^ s+ ^- £ "^\\A q p(T)\\ L P dT 

q&L q& Jo 

re[0,t] 

< CbfllpolU- + SU P / 2*'- e «M||A g /(r)||j*drY 

Proof: The proof is similar to that of Proposition 4.9 in [11] (see also Theorem 3.12 in 
[9]). First, we localize equation (|26p in the Fourier space by mean of the operator 
A q . We get 

9 t A g p + Viu • VAyp - vAA qP = A q f + 
with F q = Fl + F% + F 9 3 + F g 4 anc@ 

F l-= Yl [S g i-iu,A q ]-VA q ,p, F g 2 := ^ (5 9 _! - S g >-i)u • VA q A q ,p, 

\q'-g\<4 W-q\<l 

Fl := -A,( J] ^^^A^'), F 4 := - £ ^A^A^ £ Ag^iA 

k'^|<4 g'>g-3 ^ |q|<i ' 

Multiply both sides by | Aqp^ 1 sgn (A q p), integrate over and apply Holder 
inequality. Owing to div5 g _iu = 0, we get: 

-i\\ A qp\\ P Lp- u I ^AgplAgplP- 1 sgn^Agp) dx < (llAjlliP + ll^lli^llA^H^ 1 . 
pat j r n \ j 



2 with the summation convention over repeated indices i 
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Now, according to Lemma A. 5 in [7], there exists some constant n depending only 
on N and on Supp <p such that 

- f AA q p\A q p\P- 1 S ga(A q p)dx>K2 2 i\\A q p\\ p LP . 
Jr n ' ' 

Therefore, we end up with 

~\\\p\\l P + K u2^\\A q p\\ p LP < (\\A q f\\ LP + £ \\F*\\J) \\A q p\\%\ (29) 

i=l 

In the following calculations, assume that 2 < p < oo. Using a standard commuta- 
tion estimate for bounding F q (see e. g. [5]) and the definition of operators A q and 
S q yields 

II^IIlp< E \NS q ^u\\ L ^\\A ql p\\ LP , \\F^\\ LP < £ 2 q \\A q u\\ L ~\\A q p\\ LP , 

\q'-q\<4 W-q\<l 

\\F q 3 \\ LP < E ^HA^pll^llVlli*. E ^<^\\A q/p \\ LP \\A ql u\\ LP , 

q'<q+2 9^3-3 

with Aq := J2\ a \<N A<j' f° r some large enough positive integer No. 
Note that 

||V5,/_i«|| LO c < C E 2 9 " {1+ p ) ||A (? »u|| L p, 

q"<q'-N 

„N „ - „/ N 

||A g «|| £ =o < C2 q v \\A q u\\ LP , \\A q ,p\\ L ~> < C2 q v \\A ql p\\ L P, 



so that plugging the above inequalities in (|29j) . we get 
1 d 



^^||A 9 p||^ + ^||V||^<(||A g /|| LP 

+C E 2^'^e' q , \\A q/p \\ LP + C E 2^'-)( 1+ f )<||A^|| LP ) HA^- 1 
q'>q—i q'<q 

with e r (t) := / V 2 r ' {1+ f ) \\ A r >u\\ LP dr. 
Jo r ,< r 

Let A > be a large enough positive parameter (to be fixed hereafter). We set 

p X q (t) := 2«-e-^W||A g p(t)|| L p and / g A (i) := W> e ->*&) \\A q f(t)\\ LP . 
Obviously, the above inequality rewrites 

~(p x g r + xs' q ( P x q r + K V 2^y < (pj)f- 1 

x (/* + C2-e~^ ( E 2^'^f ^' qf \\A ql p\\ LP + E 2^'^( 1+ f )<||A^|| XP )) 



q'>q— 4 <j'<<? 
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so that, performing a time integration, we eventually get: 

p x q (t) + K v2 2 « f p x q (T)dr + \ [ t e' q (r)p x (T)dT<p x (0)+ f f X (r)dr 
Jo Jo Jo 

+C Y, 2^^ 1+ f +s ^ f e' ql {r)e x W-^p x {r)dT 

+c J- aV-tfd+f-) [\' q ( T )e x ^'^p x Ar)dT. 

Using the decomposition e' q (r) = s'/(t) + (e' q — s'i){t) in the last term and the 
fact that sequence (e n )nez is nonnegative and nondecreasing, we gather that for all 

p x (t)+ K v2 2 « f p x {T)dT + \ [ t e' q (r)p x (T)dT<p x (0)+ f f x (r)dr 
Jo Jo Jo 

+C J2 2^')( 1+ f+ s ) e ^V(*)-^W) f e' ql (r)p x (r)dT 
q'>q 

+C £ 2 (*'-*)(l+£- S ) fe' ql {r)p x (r)dr + ^Y. 2{q " q){1+ ^ ^ 

q'<q J( > ' A q ><q re[0,t] 

Suppose now that the following condition is satisfied: 

sup2 9 ( 1+ i~) ||A 9 u|| L i f£/ p\ < elog2 for some e such that \e<-( — +l+s). (30) 
qez ' 2 p 

This ensures that for all q' > g, we have 

2 ( 9 -,')(l+f +») e A( V (t)-e,(t)) < 2 ^ 1+ f +s \ 

Taking the supremum with respect to q in the equation preceding (|30|) , we thus get 
sup (p x (T) + K »2 2 * [ T p x (T')dr' + X [ T e' q (r')p x (T')dr') < sup^(0) 

q£Z \ Jo JO J q& 

re[0,i] 



+ sup / f x {r)dT + Csup / e '( r )pJ(r)dr+ y sup pj(r). 
gez j o qez Jo A <?ez 



rG[0,t] 

Therefore, 



sup p x (r) + usup2 2q / p x (r)dT + Xsup e' (r)p X (r) dr < 3||/Oo|U s 

qeZ 9 GZ Jo gGZJo p '°° 

+3sup / /g(r) fir + 3C sup / e'/(r)pg/(r) dr + — sup p g (r). 

q&L JO qeZ JO A q& 



re[0,i] 



re[0,t] 
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In order to conclude, it is only a matter of choosing A = 6C. We get 

sup P x {t) + vsu V 2 2q f Pg{r)dr < C (\\po\\bs + sup / fH^dr). 

q& qel JO V p '°° q&Jo J 



re[0,t] 

This is exactly what we wanted. 



In the case p < 2, the above bound for F q turns out to be wrong. It may be replaced 



,„ F* 

however by the following inequality 



\\Fl\\ L p<C Y, 2 9(1 V)||A g ,p|| LP ,||A (? m|| L P. 

n' ( — — — ~\ • ( — 'V1 

Hence, knowing that ||A 9 /p|| LP / < C2 V k p V \\A q ,p\\ L p, the term 2 W 9 n + p ; has 

to replaced by 2^ q q ^ 1+ V^ in all the summations over indices (q,q') such that 
q' > q — 4. Again, this leads to the desired inequality. ■ 

Let us now extend the previous estimate to the nonstationary Stokes system ([2 



Proposition 7 Let u be as in Proposition [6] and v satisfy the nonstationary Stokes sys- 
tem (J2ZJ. Then 



sup 2 qs - £ ^\\A q v(T)\\ LP + vsup f 2 q ^- £ ^\\A q v{T)\\ LP dT 

q& qeZ JO 

, s +sup f 2 qs - e ^\\A q Vf{T)\\ LP dT\ 

q&Jo J 



re[o,t] 

<C ( |b 



where V stands for the Leray projector over solenoidal vector fields. 

Proof: Applying operator V to the identity 

d t A q v + S q -iu ■ VA q v - vAA q v = A q f + F q , 

and using that VA q v = A q v, we get 

d t A q v + S q -m ■ VA q v - uAA q v = VA q f + VF q + [S^-iu, V] ■ VA q v. 

On the one hand, because TF q is supported in some annulus 2 9 C(0, r%, r-z), there 
exists a constant C > such that 

\\VF q \\ LP <C\\F q \\ LP . 

On the other hand, standard commutator estimates (see e.g. [5]) ensure that the 
new term [S q -iu, V] • VA q v satisfies the same inequalities as F q . Arguing as for 
system ([26]) . it is then easy to complete the proof. ■ 
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4.2 Proof of Theorem H 

Remind that the system satisfied by the difference between the two solutions reads 
d t S9 + div ( Ul S9) = - div (9 2 Su), 

dtdu + div (iti ® Su) — vASu + V(5Q = — div (Su (g) U2) + (59 eAr. 

We aim at proving that (SB, Su) = 0. To achieve it, we shall apply proposition [6] or [7] to 
the two equations of the above system with s = —2 + y — r\ (for some positive r\ such 
that r\ < — 1 + iVmin(|, l), which is consistent with the assumption p < 2N). 

So we first have to justify that S9 and Su belong to L^?(i?p j0 o p ) for all e E]0, 1[. 
Combining interpolation with the assumptions on the solutions (ui,0i) and (u2,9 2 ), we 

see that Ui belongs to every space L r T (Bp t00 r ) with 1 < r < 00. Because 

d t S9 = - div (6 2 u 2 - 6 lUl ), 59(0) = 0, 

Proposition H] and Holder inequality enable us to get S9 € C([0,T]; B Pi00 p ) for all 
e E]0, 1[. Plugging this new information in the equation for Su and using again Proposition 

BJ it is then easy to justify that Su is also in L^ > (B Pi00 p ) . 

One can now tackle the proof of uniqueness. Assume that the constant c has been 
chosen so small as condition ([28]) to be satisfied by the vector field u\. Denoting 

e q (t) = cY,2 q ' {1+ ^\\\>MLUL?), 

q'<q 

»(t):= sup 2-^- 2+ T-^-^\\A g S9(r)\\ LP , 

re[0,t] 

qez 

ft 

SU(t):= sup 2- q{ - 2+ T-^- £ "^\\A g Su(T)\\ LP + usnp / 2 qi f \\A q Su\\ LP dr, 

re[o,t] ?ez Jo 

q& 

we thus get according to Propositions [6] and [3 

SB(t) <Csup [ l 2 q( - 2+ f- v) - £ ^ T) \\A q div (6 2 Sa)\\ LP dr, 
q Jo 

SU(t)<Csup[ 2 q{ f ^ 2 ^- e ^ T \\\A q div (Su ® u 2 )\\L P + \\A q S9\\ LP ) dr. 
q Jo 

Let us admit that the nonlinear terms may be bounded as follows (see the proof in the 
appendix) : 

sup / 2 qi - 2+ f~ v) - eq(T) \\A q div(6 2 Su)\\ LP dT 
q Jo 

<C\\9 2 \\ sup / 2 q{ f~ v) ~ £q(T) \\A q Su\\ LP dT, 

L^(B p>x p ) q Jo 
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sup [ t 2 qi - 2+ f- r >')- e ^\\A q div(du®u 2 )\\ L pdT 



q Jo 



<C\\u 2 \\ ,n_, sup / 2 q{ f- T > ) - £ '> {T) \\AqSu\\ L pdT. 



LfiB^') q Jo 

We eventually get 

sd(t) < Cv- 1 \\e 2 \\ ._ 1+ k su(t), 

SJ(t) < cfu^l^W K_ 1 SJ{t)+ f (5B(r)dr), 
v Lf{B p p x ) Jo ' 

whence, if c has been chosen small enough in (JH), 

3U{t)<Cu- 1 [ \\e 2 \\ jv_! SJ(r)dT. 

Jo L?(B P % ) 

Gronwall Lemma thus ensures that SJ = on [0, T]. Of course, this also entails that 
dS = on [0,T]. m 

4.3 The limit case p = 2N 

Carrying out the method which has been used in the previous section to the limit case 
p = 2N seems hopeless. Indeed, we would have to deal with the product 8 2 5u while 

8 2 £ £t(-^27V oo ) an d &u £ ^tC^Iatoo) f° r some £ > 0) which does not make sense since 
the sum of indices of regularity is negative. 

In order to meet the index p = 2N, one may resort to Besov spaces with third index 1 
so as to have a velocity field in L^(Lip). Thus no losing a priori estimate is needed. On 
the other hand, due to the weak regularity assumptions, we shall be in the limit case for 
product laws in Proposition 0] so that a logarithmic interpolation argument (similar to 
that which has been used in [8] and |10| ) will be required. Let us state the result. 

Theorem 5 Let (#i,ui,VIIi) and {9 2 ,u 2 ,VH 2 ) satisfy ([1]) with the same data. Assume 
that 

9i E L§?(B~j, J and m G Lf[B^) n L^(B^ N1 ) for i = l,2. 
Then (^^VrLj = {9 2 ,u 2 ,VU 2 ) on [0,T]. 

Proof: We omit the proof of the fact that (<S9, 5u, VST) belongs to the space 

G T := Lf(BjJ x (C([0, T}; BJJ n 4( B fiv,j) x (^( Voo)) • 
In the following computations, the space Gt will be endowed with the norm 

||(0,u,vn)|| GT d = \\e\\ 3 +||wi| 3 +|HU * +l|vni|_ 3 . 

T ^AT.ooJ T ( a 2N,oo> T \ 2iV,oo / L, T\ ts 2N,ao> 
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A priori estimates for the transport equation (see e.g. the limit case in Proposition 
4.7 of [11]) and inequality (jlip guarantee that for all t < T 

m L~(B- 3 * <^e X p(c||V Ml || , ) [\\H Bi INI J cJr. (31) 
Next, we have according to Proposition 3.2 of [2] and Inequality (|16p . 
INI 3 +H|«5u|L i +||Vdn|L 3 

(^JV.oJ L t( B 2N,oo> L t( B 2N,oo) 

C|l«l|l x 3 

<Ce Ll t {B 2N,i) \\S9e N - div(6u®u 2 )\\ dr. 

Because div&i = 0, we have for 1 < i < iV, 

( div (& u 2 )) 4 = dj (f^nl + uj)) + 2^ 

So using Proposition [3l we find that 

II div (5a®u 2 ) II ._8 <C||&t||._3 1 1 1 1 . a . 

ft ? D p>^ 

a 2N,oo -°2JV,oo - D 2JV,1 

Now, taking advantage of Proposition 1.8 in [SJ, one may write 



x ._i + 

/,| C I 1 / L ti B 2N,l> L t( B 2N,l) 

i <C\\6u\\_ i log eH ——. — 

L t ( B 2N,oo> 

Let us introduce the notation 

V(t) = t[\\Ux\\ ,_1 +\\U 2 \\ ._! )+||«l|| 1 .3 +||«2|| 1 .f , 

L T( B 2N,l) L T( B 2N.l) L t( B 2N,l) L t( B 2N,l> 

W{t) = ^h-lB-* ) + M I^ ) + l|VOT|l &fA-* )■ 

L t ( B 2N,oo> L t( B 2N.oo> L t\ B 2N,oo> 



Because 



lj t\ t >2N,l> ^t\ ti 2N,ll 



and the map x i — ► xln(e + ^) (for fixed y > 0) is nondecreasing over M + , we end 
up with 

C||Ul|| 3 / -rr/ \ \ 

W(t)<Ce / (||n 2 ( S )||.3 -f||0 2 (s)||. _i ) WW log e + -f( ) ds. 



Applying Osgood lemma (see e.g. [5]) thus yields W = on [0, T] whence also 
8) = according to inequality ([3T|) . This completes the proof. ■ 
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5 Appendix 



This appendix is devoted to proving the estimates for the convection terms that we used 
in section B~2l 



Lemma 4 Let (a q ) q£ z be a sequence of nonnegative functions over [0, T]. Let si, S2, p 
satisfy 

N N ( 2 \ 

1 < P < oo, V 1 > s±, — > S2 and si + S2 > iVmax 0, 1 . 

p p V p J 

Assume that for all q' > q and t E [0, T], we have 

0<a q ,(t)-a q (t) < U Sl + s 2 + Nmm(0,l-~)\q ~q)- (32) 

Then for all r € [l,oo], there exists a constant C depending only on s±, S2, N and p 
such that for all function b and solenoidal vector field a over M. N , the following estimate 
holds true for all t € [0, T] ; 



sup / 2 

q& JO 



A g div(a6)||iP dr < C\\b\\r r( ^i sup 

tK p '°°' q&Z 



2 ^-^\\A q a\\ LP 



Proof: The proof relies on Bony's decomposition. Knowing that div6 = 0, we have 
(with the usual summation convention over repeated indices): 

A q div (ab) = A q {T djh <j) + A q (t a ,d 3 b) + A q djR{a^ b). (33) 

By virtue of ([7]), one may write 

A q (f a . b ai) = \iStf-idjbAjai). 

\q'-q\<4 

For the sake of simplicity, let us proceed as if AqiTg.^a^ = S q -\djb A q a^ (having 
([32]) justifies this approximation). Using the definition of S q -i, we get 



2 q ( Sl - 1 -f)\\S q _ 1 d j bA q a j \\ LP < || A 



q'<q-2 



''''^>||A,V6||^2^<*^ 



In consequence, we have for all < t < T, 

J' ^••+«- 1 -T>-^>\\S t _ 1 a j b^a>\\ L ,dT 



2<? S2-« 9 || A II p 



^^||A (? ,V6|| L , (L o 0) 2^')( Sl - 1 -f). 



q'<q-2 



Since s\ — 1 — ^ < 0, we thus get 



sup / 2 



1 q(s 1 +s 2 -l-^)-a q {T) I 



A g T 9 4 hOp\\LP dr 



<C||V6||_ . Sl -!-K sup 

^( S °o,oo p ) qez 



2^- a "\\A q a\\ LP 



(34) 
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Likewise, in order to bound the second term of ([33]) . one may proceed as if 
Now, for all < r < T, we have 

2^ sl+S2 - 1 -f ) - a " iT) \\s q ^A q d j b\\ LP 

<2^-^\\A q Vb\\ LP 2- q ' { ^- S2) - a « {T) 

q'<q-2 

Knowing that a„ > a q > for q > q', we deduce that 



> <?(si+S2~l-^)-a 9 (r) 



\\S q ^a? A q dM LP dr 



q'<q-2 



'||Ag/0 



2 (9'-9)(f-«). 



Combining Bernstein inequality and the fact that S2 < ^, we thus conclude that 



sup /V*^- 1 - 
ogz Jo 



: )-a,(r) 



\\S q ^a 3 A q d 3 b\\ LP dr 

< C\\Vb\\r r(BS1 -i, sup 2^-^\\A q a\\ LP 



(35) 



In order to treat the remainder term, we shall consider the cases p > 2 and p < 2. 
Let us start with the case p > 2 which is slightly easier. We have 

A q djR(a j ,b)= Y djA q (A q/ a j A q >b) with A g / = A ? /_i + A g / + A q > + i. 

q'>q-3 

Because T^djAg^Aq/a 3 'A q >b)) is supported in a ball of size 2 q , Bernstein inequality 
ensures that 



->gOl+s 2 -l--^)-a q (T) 



WA^R^MW^C Y, 2^ + ^~ a ^\\A q ,a\\ LP \\A q ,b\\ LV , 



q'>q-3 



whence 

f 2 9(si+S2 - 1 -f ||A g %i?(^, 6)||iP dr 

J o 

<C V / 2 g ' S2 ^'Wj|A^a|| L P (2 9 ' Si ||A 9 ,6|| Z p) 2K'^)M 2 (9-?' 

Thanks to assumption (|32|) . we have 

(<V - a ? )(r) + {q- q'){si + s 2 ) < (q - q) (-^y- ^) for a11 4 ^ <?• 
As si + S2 > 0, we end up for all q € Z with 



)(S1+S2) 



,g(si+s 2 -l-— )-a 9 (r) 







tK p ' 00 > q & 



which, together with 



and ([3"5"j) . completes the proof in the case p > 2. 
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If p < 2, one may write (use Bernstein inequality) 

2 q{s ^ 1 -^- a ^\\A q d j R(a j ,b)\\ LP <C 2 9(si+S2 ^-f ) - a9(T) ||A,,a|| LP H|A (? ,6|| L p 

q'>q-3 
q'( — - N ) ■ 

and use that ||A g /a|| iP / < C2 ^p 7' \\/\ q ia\\LP- Therefore, one may go along the 
lines of the case p > 2. It is only a matter of changing the term 2^ q ~ q )( s i+ s 2) m t Q 

2 (,-g')( il + a2 + £-£)_ " m 
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